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Abstract 

In this work we derive an approximated no-arbitrage market valuation formula 
for Constant Maturity Credit Default Swaps (CMCDS). We move from the CDS 
options market model in Brigo (2004), and derive a formula for CMCDS that 
is the analogous of the formula for constant maturity swaps in the default free 
swap market under the LIBOR market model. A "convexity adjustment" -like 
correction is present in the related formula. Without such correction, or with zero 
correlations, the formula returns an obvious deterministic-credit-spread expression 
for the CMCDS price. To obtain the result we derive a joint dynamics of forward 
CDS rates under a single pricing measure, as in Brigo (2004). Numerical examples 
of the "convexity adjustment" impact complete the paper. 

Keywords: CDS Options, CDS Options Market Model, Constant Maturity CDS, 
Convexity Adjustment, Participation Rate, CDS rates volatility, CDS rates correlation. 
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1 Introduction 

Constant Maturity Credit Default Swaps (CMCDS) are receiving increasing attention in 
the financial community. In this work, we aim at deriving an approximated no-arbitrage 
market valuation formula for CMCDS. We move from the CDS options market model 
in Brigo (2004), and derive a formula for CMCDS that is the analogous of the formula 
for constant maturity swaps in the default free swap market under the LIBOR market 
model. 

In Brigo (2004) the focus is somehow different, since in that paper we derive CDS 
option prices for CDS payoffs given in the market and for some postponed approximated 
CDS payoffs. We do so by means of a rigorous change of numeraire technique following 
Jamshidian (2002) and based on a result of Jeanblanc and Rutkowski (2000). We also 
establish a link with callable defaultable floaters. 

This paper starts by recalling again one alternative expression for CDS payoffs, stem- 
ming from different conventions on the premium flows and on the protection leg. We 
briefly introduce CDS forward rates, postponing the detailed definitions. We define 
CMCDS and give immediately the main result of the paper, the approximated pric- 
ing formula, in terms of CDS forward rates and of their volatilities and correlations. 
We point out some analogies with constant maturity swaps, showing that a "convexity 
adjustment "-like correction is present. Without such correction, or with zero corre- 
lations, the formula returns an obvious deterministic-credit-spread expression for the 
CMCDS price. 

Once the main result has been described, we move to introduce the formal apparatus 
that allows to prove it. We introduce rigorously CDS forward rates. This leads to an 
investigation on the possibility to express such rates in terms of some basic one-period 
rates and to a discussion on a possible analogy with the LIBOR and swap default-free 
models. We discuss the change of numeraire approach to deriving a joint dynamics of 
forward CDS rates under a single pricing measure, derivation that is only hinted at 
in Brigo (2004) for deriving the Black-like formula for CDS options. Through a drift- 
freezing approximation we then prove the formula for CMCDS pricing and give some 
numerical examples highlighting the role of the participation rate and of the convexity 
adjustment. 

Remark 1.1. (How to use this paper) The CDS Options market model is the same 
as in Brigo (2004) '^^^ based on a complex change of numeraire involving two families 
of forward CDS rates. This apparatus is reported in Section^hut there is no need to go 
through it if one does not plan to test the formula against Monte Carlo simulation. As 
far as the formula derivation is concerned, the reader may skip Section^ and go directly 
to the much simpler approximated model in Section 
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2 Credit Default Swaps (CDS) and Constant Matu- 
rity CDS 

We recall briefly some basic definitions for CDS's. 

Definition 2.1. (Credit Default Swap). A CDS is a contract where the "protection 
buyer" "A" pays rates R ("R" stands for Rates) at times Ta+i, ... ,31, (the "premium 
leg") in exchange for a single protection payment Lgd (loss given default, the protection 
leg). "A" receives the protection leg by the "protection seller" "B" at the default time r 
of a reference entity "C", provided that Ta < t < T^. The rates R paid by "A" stop in 
case of default. We thus have 



Protection Seller — > protection Lgd at default tc of "C" if in [T^, T^] Protection Buyer 
"B" ^ rate R at r„+i, . . . , Tb or until default tc ^ "A " 

This is called a "running CDS" (RCDS) discounted payoff. 

We explicitly point out that we are assuming the offered protection amount Lgd to 
be deterministic. Typically Lgd = 1 — Reg, where the recovery rate Reg is assumed to 
be deterministic and the notional is set to one. 

Sometimes a slightly different payoff is considered for RCDS contracts, fnstead of 
considering the exact default time r, the protection payment Lgd is postponed to the first 
time Tj following default. If the grid is three or six months spaced, this postponement 
consists in a few months at worst. 

We term "Postponed payments Running CDS" (PRCDS) the CDS payoff under the 
postponed formulation. The advantage of the postponed protection payment is that no 
accrued-interest term is necessary, and also that all payments occur at the canonical grid 
of the Tj's. The postponed payout is better for deriving market models of CDS rates 
dynamics, as we shall see shortly. It is also fundamental in establishing the "defaultable 
floater analogy" as we have seen in Brigo (2004). 

We denote by CDSa_b(t, i?, Lgd) the price at time t of the above standard running 
CDS flows to the protection seller "B" . We add the preflx "PR" to denote the analogous 
price for the postponed payoff. The pricing formulas for these payoffs depend on the 
assumptions on interest-rate dynamics and on the default time r. In general the CDS 
forward rate i?a^b(0) for protection in T^, at initial time is obtained (in the postponed 
case) by solving the equation PRCDSa,6(0, i?a,fe(0), Lgd) = 0. We will give details on 
this later: now we only say that the market provides quotes for i?o,6(0)'s for increasing 
maturities T5 (notice that Tq = 0, so that only "spot" CDS rates are quoted in the 
market). The i?a,6(0) rate makes the CDS contract fair at the valuation time. A special 
role in our work will be covered by one-period rates i?i(0) = i?i_i,i(0) (protection in 
[Tj_i, Tj]). These rates may seem artiflcial but they can be easily computed from quoted 
spot rates. An important role is also assigned to the corporate zero-coupon bond price 
P(0, T), which is the price at time of one unit of currency made available by name "C" 
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at maturity T if no default occurs, and with no recovery in case of early default. The 
corresponding default-frcc zero coupon bond is denoted by P{0,T). 

We are now ready to introduce CMCDS's. 

In a CMCDS with first reset in Ta and with final maturity Th, protection Lgd on a 
reference credit "C" against default in [Ta,Tfe] is given from a protection seller "B" to a 
protection buyer "A" . However, in exchange for this protection, a "constant maturity" 
CDS rate is paid. 

We know by definition that the fair rate to be paid at Tj for protection against default 
in [T'i_i,T'j] would be Ri. This leads us to the following 

Remark 2.2. (A "floating-rate" CDS). A contract that protects in Ta,Ti, can be 
in principle decomposed into a stream of contracts, each single contract protecting in 
[Tj-i, Tj] , for j — a + 1, . . . ,b, say with protection payment Lgd postponed to Tj if default 
occurs in [Tj_i ,Tj]. In each single period, the rate Rj (Tj^i) paid at Tj makes the exchange 
fair, so that in total a contract offering protection Lgd on a reference credit "C" in 
[Ta, Tb] in exchange for payment of rates Ra+i{Ta), . . . , RjiTj^i), . . . , Rb{Tb-i) at times 
Ta+i, . . . ,Tj, . . . ,Tb is fair, i.e. has zero initial present value. This product can be seen 
as a sort of floating rate CDS. 

However, in CMCDS's the rate that is paid at each period for protection is not the 
related one-period CDS rate, as would be natural from the above remark, but a longer 
period CDS rate. Consider indeed the following 

Definition 2.3. (Constant Maturity CDS). Consider a contract protecting in [Ta,Tb\ 
against default of a reference credit "C". If default occurs in [Ta,Tb], a protection payment 

Lgd is made from the protection seller "B" to the protection buyer "A" at the first Tj 
following the default time. This is called "protection leg". In exchange for this protection 
'A" pays to "B" at each Tj before default a "c + 1-long" (constant maturity) CDS rate 
Rj-ij^c{Tj-i) (times a year fraction Uj = Tj — Tj_i), with "c" an integer larger than 
zero. Notice that for c = we would obtain the fair "floating rate" CDS above, whose 
initial value would be zero. 

Given that c > in our definition, the value of the contract will be nonzero in general, 
so that we have to find this vahic at the initial time if we arc to price this kind of 
transaction. We face this task by resorting to the market model derived below. 

The value of the CMCDS to "B" is the value of the premium leg minus the value of 
the protection leg. The protection leg valuation is trivial, since this is the same leg as in 
a standard forward start [T^, T^ CDS. As such, it is for example equal to 

h b 
j=a+l j=a+l 

This value has to be subtracted to the premium leg. The non-trivial part is indeed 
computing the premium leg value at initial time 0. Notice that the final formula can be 
implemented easily on a spreadsheet, requiring no numerical apparatus. 
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Proposition 2.4. (Main Result: An approximated formula for CMCDS) Con- 
sider the Constant Maturity CDS defined in Definition \2.3[ The present value at initial 
time of the CMCDS to the protection seller "B" is 

CDSa.^,AO, Lao) = f: «,^(0, T,) (1) 



Ri{0) exp 



J=a+1 V i=J 

cTkRkiO) 

■' P. / 



k' 

RkiO) + Lao/ak 




where Rk{0) are the one-period CDS forward rates for protection in [Tk^i,Tk]. These 
CDS forward rates can be computed from quoted spot CDS rates -Ro,fc(0) and corporate 
zero coupon bonds P{0,Tk) via 



P ^0.^(0) Eti (^hPjO, n) - Ro,k-i{0) Etl (^hPjO, T,) 

while Rk{0) are approximations of the Rk{0) (equal in case of independence of interest 
rates and credit spreads ) in terms of corporate P and default free P zero coupon bonds 
given by 

RkiO) . ^.(0) ^ ,,,mT^P(0,T.)/P(0,T._0 - P(0,T.) 



akPiO,Tk) 

and where: ak is the volatility of Rkit), assumed constant (we deal with the time-varying 
volatility in the proof below); 

Pi j is the instantaneous correlation between Ri and Rj; 

One-period forward CDS rates volatilities ak can in principle be stripped from longer 
period CDS volatilities, similarly to how forward LIBOR rates volatilities can be stripped 
from swaptions volatilities in the LIBOR model. This stripping is made possible from an 
approximated volatility formula based on drift freezing (formula (6.58) for the LIBOR 
case in Brigo and Mercurio (2001)). Cascade methods are also available for this (as in 
Brigo and Morini (2004)), although for the time being the only available CDS options all 
have short maturities and the lack of a liquid market discourages this kind of approach. 
For the time being the above formula can be employed with stylized values of volatilities 
to have an idea of the impact of the "convexity adjustments" . Finally, one may consider 
using historical volatilities and correlations in the formula as first guesses. 

As a further remark we notice that, if not for the exponential term (which vanishes 
for example when p's are set to zero) this expression would be, not surprisingly. 



CDScMa,bAO, Lgd- p = 0) = a,P(0, T,)(P,_i,,+e(0) - Rj^iM) (2) 

j=a+l 
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Notation 


Description 


Ta, {Ta+1, ■ ■ ■ , Tb-l), Tb 

ai 
L{S,T) 


Default time of the reference entity "C" 
initial and final dates in the protection schedule of the CDS and CMCDS 
First of the Tj's following t and i, respectively 
year fraction between Ti_i and 
LIBOR rate at time S for maturity T 


Ra,b 
Ri 

Ri-2,i 

Rec 
Lgd = 1 — Reg 

nPRCDSa,fc(t) 

PRCDS„,b(t,i?, Lgd) 
CDScMq (, ,,(^7 Lgd) 


Rate in the premium leg of a CDS, paid by "A", the protection buyer 
CDS Rate to be paid by "A", the protection buyer at for protection in \Ti^i^Ti\ 
CDS Rate to be paid by "A" at Ti_i and Ti for protection in [Ti_2, 7i] 
Recovery fraction on a unit notional 
Protection payment against a Loss (given default of "C" in \Ta,Tb]) 
Discounted payoff of a postponed running CDS to "B" , the protection seller 
Price of a running postponed CDS to "B", protecting against default of "C" in [T^jTb] 

Price at time t of a CMCDS to "B" , protecting against default 
of "C" in [Ta,Tb] in exchange for a periodic constant maturity "c+ l"-long CDS rate 


l{r>T} 
l{r<T} 
B{t) 

D{t,T)^B{t)lB{T) 
P{t,T) 

l[r>t}P{t,T) 
Gt 

DC{-) 


Survival indicator, is one if default occurs after T and zero otherwise 
Default indicator, is one if default occurs before or at T, and zero otherwise 
Bank account numeraire of the risk neutral measure at time t 
Stochastic discount factor at time t for maturity T 
Zero coupon bond at time t for maturity T 
Defaultable Zero coupon bond at time t for maturity T 
Defaultable "Present value per basis point" numeraire and associated measure 

Default free market information up to time t 
Default free market information plus explicit monitoring of default up to time t 
DC{Xt) is the row vector v in dXt = {...)dt + v dWt for diffusion processes X 
with W vector Brownian motion common to all relevant diffusion processes 



Table 1: Main notation in the paper. 



The exponential term can be considered indeed to be a sort of "convexity adjustment" 
similar in spirit to the convexity adjustment needed to value constant maturity swaps 
with the LIBOR model in the default-free market. 

Finally, this formula should be tested against prices obtained via Monte Carlo simu- 
lation of the dynamics DJ]) before being employed massively. One should make sure that 
for the order of magnitude of volatilities, correlations and initial CDS rates present in 
the market at a given time the freezing approximation works well. We plan to analyze 
this approximation against Monte Carlo simulation in further work. This future work is 
the reason why we derive the exact rates dynamics below. 
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3 CDS Option Market Model Dynamics 

In the remaining part of the paper we build the apparatus allowing us to prove the main 
result rigorously. We specify the probabilistic framework in the following remark. 

Remctrk 3.1. (Probabilistic Framework: r as first jump of a Cox process) 

Here we place ourselves in a probability space {Q, Q, Q) where the default time random 
variable r will be defined. The probability measure Q is the risk neutral probability. 
This space is endowed with a filtration {J^t)t, typically representing the basic filtration 
without default, i.e. the "information flow" of interest rates, intensities and possibly 
other default-free market quantities. Obviously J^t Q G for all t. 

We consider a non-negative, {J-t)t progressively measurable process A with integrable 
sample paths in (Q, Q, Q). 

The space {Q, Q, Q) is assumed to be sufficiently rich to support a random variable U 
uniformly distributed on [0, 1] and independent of {J-t)t- The random default time r can 
then be defined as 

T := inf |i > : exp Kds^ < 

With this definition A is indeed the Tt stochastic intensity of the default time t, in 
that Q(t > t\J^t) = exp(- /J Xsds). 

We consider also the filtration Qt — V (7({t < s}, s < t) and assume Gt Q G- 
The second sigma field cr({T < s},s < t) contributing to Qt represents the information 

on whether default occurred before t, and if so when exactly. Since this information is 
available to us when we price, we need to condition on Qt rather than on Tt alone. 

For more details on the canonical construction of a default time with a given hazard 
rate see e.g. Bielecki and Rutkowski (2001), p. 226. 

The remark above amounts to saying that default r is modeled as the first jump time 
of a Cox process with the given intensity process. We will not model the intensity directly 
but rather some market quantities embedding the impact of the relevant intensity model 
that is consistent with them. An explicit tractable stochastic intensity/ interest-rate 
model with automatic analytical and separable calibration to interest rate derivatives 
and CDS's is given for example in Brigo and Alfonsi (2003), where an analytical formula 
for CDS options based on Jamshidian's decomposition is also presented. 

Formally, we may write the RCDS discounted value at time t as 

nRCDStt,b(i) := DiscountedPremiumLeg — DiscountedProtectionLeg (3) 

b 

= D{t,T){T - Tf3^r)-l)R'i-{Ta<T<n} + D{t,Ti)aiRl{r>Ti} - l{Ta<r<n}D{t, t) Lgd 

i=a+l 

where t e [T^(t)_i, T^(t)), i.e. T^(t) is the first date among the Tj's that follows t, and 
where ctj is the year fraction between Tj_i and Tj. The stochastic discount factor at time 
t for maturity T is denoted by D{t, T) — B{t)/B{T), where B{t) denotes the risk-neutral 
measure bank-account numeraire. 
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Under the postponed formulation, where the protection payment is moved from r to 
T^(t-), the CDS discounted payoff can be written as 

b b 

npRCDSa,fe(t) := ^ D{t,Ti)aiRl{r>Ti} - ^ l{T,_i<T<r,}-D(t, Tj) Lgd, (4) 

i=a+l i=a+l 

which we term "Postponed payments Running CDS" (PRCDS) discounted payoff. Com- 
pare with the earher discounted payout ([3]) where the protection payment occurs exactly 
at r: The advantage of the postponed protection payment is that no accrued-interest 
term in (r — T^(t-)_i) is necessary, and also that all payments occur at the canonical grid 
of the Tj's. 

In general, we can compute the CDS price according to risk-neutral valuation (see 
for example Bielecki and Rutkowski (2001) for the most general result of this kind): 

CDSa,b{t, R, Lgd) = E {UnCUSa,b{t)\Gt} (5) 

where we recall that Qt = ■^t^ ^({^ < u},u < t). In the Cox process setting default 
is unpredictable and this is why observation of !Ft alone does not imply observation of 
the default time, contrary to standard structural (Merton, Black and Cox, etc) models 
where instead J^t = Qt- At times we denote by and Qt the expectation and probability 
conditional on the default-free sigma field J^f The above expected value can also be 
written as 

CDS,,,(t, R, Lgd) = ^^I^Le {nHCDS,,,(t) l^i} (6) 

(see again Bielecki and Rutkowski (2001), or more in particular Jeanblanc and Rutkowski 
(2000), where the most general form of this result is reported). This second expression, 
and especially the analogous definitions with postponed payoffs, is fundamental for intro- 
ducing the market model for CDS options in a rigorous way. We explicit the postponed 
expression by substituting the payoff: 

PRCDS„,fe(t,i?) = -i^{-LGD MMT,-,<r<T,}D{t,T,)] + R J] Et[Z)(t, rOaa{.>T,}] } (7) 

i=a+l i=a+l 

Let us deal with the definition of postponed (running) CDS forward rate Ra,bit). This 
can be defined as that R that makes the PRCDS value equal to zero at time t, so that 

PRCBSa,bit, Ra,b{t), Lgu) = 

(notice that this Ra,b is the of Brigo (2004)). The idea is then solving this 
equation in Rabbit). In doing this one has to be careful. It is best to start moving 
from expression ([6]) rather than ([5]). Equate this expression to zero and derive R corre- 
spondingly. Strictly speaking, the resulting R would be defined on {r > t} only, since 
elsewhere the equation is satisfied automatically thanks to the indicator in front of the 
expression, regardless of R. Since the value of R does not matter when {r < t}, the 
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equation being satisfied automatically, we need not worry about {r < t} and may define, 
in general, R =ProtectionLegValue/PremiumLeg Value, 

^ _ Loo Ela-,inDit,m{T.^,<r<T.}\J't] 

j:la^iCyMr>t\J^t)P{t,T,) 



where 



Hr>t} 



P{t,T) ■.= ¥.[D{t,T)l{r>T}\Gt] = l{r>t} E[D{t,T)l{r>T}\J't]/Q{r > t\J^t 



is the price at time t of a defaultable zero-coupon bond maturing at time T. The corre- 
sponding default free zero coupon bond is denoted by P(t,T). This approach amounts 
to equating to zero only the expected value part in and in a sense is a way of privi- 
leging "partial information" !Ft expected values to "complete information" Qt ones. Our 
R above is defined everywhere and not only conditional on r > t. The technical tool 
allowing us to do this is the above-mentioned Jeanblanc Rutkowski (2000) result, and 
this is the spirit of part of the work in Jamshidian (2002). 

A remark on how the market quotes running CDS prices is in order at this point. 
First we notice that typically the T's are three- months spaced. Usually at time t = 0, 
provided default has not yet occurred, the market sets _R to a value Rff{0) that makes 
the CDS fair at time 0, i.e. such that CDSa,fe(0, i?^|,°(0), Lgd) = 0. Actually, bid and 
ask levels are quoted for R. Typically in quoted CDS we have = and Tf, spanning 
a set of increasing maturities, even though in recent times the quoting mechanism has 
changed in some respects. Indeed, the quoting mechanism has become more similar to 
the mechanism of the futures markets. Let be the current time. Maturities Ta, . . . ,Tb 
are fixed at the original time to some values such as ly, 2y, 3y etc and then, as time 
moves for example to t = Iday, the CDS maturities are not shifted correspondingly of 
1 day as before but remain ly,2y etc from the original time 0. This means that the 
times to maturity of the quoted CDS's decrease as time passes. When the quoting time 
approaches maturity, a new set of maturities are fixed and so on. A detail concerning 
the "constant maturities" paradigm is that when the first maturity is less than one 
month away from the quoting time (say 0), the payoff two terms 

are replaced by 

Ta+lD{0, Ta+l)Rl{r>Ta+i} 

in determining the "fair" R. If we neglect this last convention, once we fix the quoting 
time (say to 0) the method to strip implied hazard functions is the same under the 
two quoting paradigms. The same happens when not neglecting the convention if we are 
exactly at one of the "0 dates", so that for example Ti — t = ly. Brigo and Alfonsi (2003) 
present a more detailed section on the "constant time-to-maturity" earlier paradigm, and 
illustrate the notion of implied deterministic intensity (hazard function) . 
We also set 

b 

j=a+l 
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the denominator in (IHl), so that Ra,b is a tradable asset (upfront CDS) divided by Cafiif). 
It follows that under the measure Q"''' having Ca,b{t) as numeraire, Ra,b follows a mar- 
tingale. For more details see Brigo (2004). We will be interested in the particular cases 
a = i — l,b = i and a = i — 2,b = i. 

In the following it will be useful to consider a running postponed CDS on a one-period 
interval, with = Tj^i and = Tj. We obtain, for the related CDS forward rate: 

J,,,, , E[D(t,r,)i|T,,,<.<T,}|^.] ^ E[D(t,r,)i{T,_,<.<T,}|^t] 

'■= L/GD — — - — — — — = L/C 



GD- 



{t > t\J^t)P{t,T,) C,_i,,(t) 



where we have set Rj := Rj-i,j- 

A last remark concerns an analogy with the default-free swap market model, where 
we have a formula linking swap rates to forward rates through a weighted average. This 
is useful since it leads to an approximated formula for swaptions in the LIBOR model, see 
for example Brigo and Mercurio (2001), Chapter 6. A similar approach can be obtained 
for CDS forward rates. It is easy to check that 

Z^i=a+l '^j-' I*' V i=a+l i=a+l 

A possible lack of analogy with the swap rates is that the w's 



ELa+1 (^hP{t, Th) 



cannot be expressed as functions of the i?i's only, unless we make some particular as- 
sumptions on the correlation between default intensities and interest rates. However, if 
we freeze the w's to time 0, which we have seen to work in the default-free LIBOR model, 
we obtain easily a useful approximated expression for Ra^b and its volatility in terms of 
RiS and their volatilities/correlations. A similar approach is pursued in Section [5] 
below, and the reader who is interested only in the derivation of the approx- 
imated formula given in the beginning may go there directly. Here instead 
we hint at deriving the real dynamics without compromises, which will be 
useful in future work for Monte Carlo tests of the numerical approximation. 

In general, when not freezing, the presence of stochastic intensities besides stochastic 
interest rates adds degrees of freedom. Now the P's (and thus the tD's) can be determined 
as functions for example of one- and two-period rates. Indeed, it is easy to show that 

r>U ^\ r>U rr x "j-l(-Ri-l(^) " -Ri-2,i(^)) P{t,Tj) -pi- Rk-i - Rk-2,k .^„x 
^{t,Ii - r'[t,Ii-i — — — — -- — , ^— — — - — II — — [iU 

ai{Ri-2,i{t) - Ri[t)) P{t,Ti) aj ^^.-^^ Rk-2,k - Rk 

To have the formula working we need to assume Ri-2.i{i) 7^ Ri{t)- We will therefore 
assume in the paper that at the initial time i?i_2,j(0) 7^ Ri{0). Under the approximated 
frozen dynamics for these two quantities we will derive below, we can see that the proba- 
bility of them to be equal at any future time t is generally 0, since this can be computed 
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as the probability of the difference of two correlated continuous random variables to be 
zero. 

We show below how this formula helps us in obtaining a market model for CDS rates. 
For the time being let us keep in mind that the exact weights w{t) in (l9l) are completely 
specified in terms of -Rj(t)'s and i?i_2,j(t)'s, so that if we include these two rates in our 
dynamics the "system" is closed in that we also know all the relevant P's. The difference 
with the LIBOR/Swap model is that here to close the system we need also two-period 
rates. 

We close this section by summarizing our notation in Table ([I]). 

4 One- and Two- Period CDS Rates joint Dynamics 
under a single pricing measure 

Let us postulate the following dynamics for one- and two- period CDS forward rates. 
Recall that Rj = Rj-ij- 

dR,{t) = a,{t)R,{t)dZi{t) 

dR,^,,{t) = u,it;R)R,.2At)dVr''\t) 

In the Brownian shocks Z and V the upper index denotes the measure (i.e. the 
measure associated with the numeraires Cj-ij,Cj^2,j in fhe above case) and the lower 
index denotes to which component of the one- and two- period rate vectors the shock 
refers. The volatilities a are deterministic, whereas the z/'s depend on the 
one-period i?'s. We assume correlations 

dZidZj = Pijdt, dVidVj = i]ijdt, dZidVj = Oijdt 

and i?i_2,i(t) G (min(i?i_i(t),[i?i_i(t) + i?i(t)]/2),max(i?i_i(t),[i?i_i(t) + i?i(t)]/2)). 
This latter condition ensures that the resulting P from formula ( 1101) be posi- 
tive and decreasing with respect to the maturity, i.e. < P{t,Ti)/ P(t,Ti^i) < 1. 
The specific definition of u ensuring this property is currently under investi- 
gation. 

We aim at finding the drift of a generic Rj under the measure associated with C*j_i^j, 
let us say for j > i. 

The change of numeraire toolkit provides the formula relating shocks under Cj_i,i 
to shocks under Cj-2j, see for example Formula (2.13) in Brigo and Mercurio (2001), 
Chapter 2. We can write 



d 



d 



Z' 



— CorrMatrix x Vector D iff usionCoefficient I In I ^ ^'•^ 1 1 dt 



Let us abbreviate "Vector Diffusion Coefficient" by "DC" 
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This is actually a sort of operator for diffusion processes that works as follows. 
DC (Xt) is the row vector v in 



dXt = {...)dt + v d 



Zt 
Vt 



for diffusion processes X with Z and V column vectors Brownian motions common to 
all relevant diffusion processes. This is to say that if for example dRi = aiRidZl, then 

DC{Ri)^[aiRi, 0, 0,..., 0]. 

Let us call Q the total correlation matrix including p, rj and 9. We have 



d 



d 



g DC In 



dt 



Now we need to compute 



DC In 



DC In 



aiP{t, Ti) 



= DC In 



Otj-l CCj T-r Rk-1 — Rk-2,k _|_ Oij Oti -pr Rk-1 — Rk-2,k 

ai Oij-i Rk-2,k ~ Rk OLi Oi'i . ^ Rk-^. k — Rh 

■' k=i+\ ' 



= DC l^ln 

=Dcf x: 1^ 



n 



Rk~i — R. 



'k-2,k 



R 



■k-2,k 



\k=i+l 



,fc=i+l 
Rk-l — Rk-2,k 



Rk 



Rj-1 - Rj-2,j 



Rj-2,j - Ri 



R 



■k-2,k 



— Rk 



+ DC In 



Rj-i - Rj 



R-j 



gDc(ln(^ 



Rk-1 ~ Rk-2,k 
Rk 



+ DC In 



'j-2,j -"-J 

Rj-1 - Rj 



Ri— 



Ri 



= J2 [DC(ln(i?,_i - Rk-2,k)) - DC(ln(i?,_2,fe - Rk))] + 



k=i+l 



+DC(ln(i?^_i - Rj)) - DC(ln(i?j_2j - Rj)) 

DC( Rk-1 — Rk~2,k) _ DC( Rk-2,k — Rk) _|_ 

Rk-1 — Rk-2; 



fe=i+l 



Rk 



Rk 



k=i+l "k-2,k 

DC( - Rj) DC( Rj-2,j - Rj) 



R,-i^R, 



Rj-2,j - Rj 



^ y ( BC{Rk-i) DC{Rk-2,k)) _ y [_DC(7?fc_2,fe) - DC(i?fe)) 

Rk-1 — Rk-2,k 



k=i+l 



k=i+l 



R, 



k-2,k 



Rk 



DC(i?,_i) - DC{Rj) DC{Rj-2,j) - BC{Rj) 



Rj-1 - Rj 



Rj-'^j ~ Rj 



It follows that 
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j^j-2,j _ _ _ Sr~^ (Pfc-l,mCfc-l-Rfc-l ~ Om,k^kRk-2,k) (^m,fc^fc-Rfc-2,fc — Pk.mCkRk) 

k=i+l ' fc=j+l 

Pj-l,m'^j-lRj-l ~ Pj,m'^jRj ^m.,jl^jRj-2,j ~ Pj,mO'jRj 



dt + '"'"^ ^"^'^ ^^''"^ ^ (it 



and 



jyj-2,j _ jyi _ _ Sr^ {dk-l,ni<^k-lRk-l — ' I niJ.-1^k:R k—2,k) jV in ,k1^kR k-2,k — k,m<^ kRk) 

-n-fe-i — J^k-2,k ^k-2,k — -n-fe 

k=i+l ' k=i+l ' 

0j-l,m'^j-lRj-l — ^j,m(^jRj_^^ _|_ Vj,m^jRj-2,j ~ ^j,m'^jRj _. TiJ 



Rj-i - Rj Rj-2,j - Rj 

Therefore, by subtracting from the first equation, taking h > i: 
~ dZl^ — dZl^ ^'^ — dZl^ — {dZl^ ^'^ — dZ!^) — 

E{Pk-l,m<^k-lRk-l — Srn,k^kRk-2,k) ST^ {dm,k^kRk~2,k — Pk,mO'kRk) _ i. 

, . , , -Kfc-l — Kk-2,k J^k-2,k — J^k 

k=i+l ' k=i+l ' 

so that we finally obtain (taking h — j) 

dRj{t) = ajRj{t){p^fdt + dZ}{t)) 
dRj^2At) = i^jRj-2AtWydt + dVj{t)), 

or, by setting 
we have 



dRj{t) = Rj{t){pidt + (JjdZ]{t)), dRj_2,j{t) = Rj_2,j{t){(t)]dt + i^jdVj{t)), 

and since p and arc completely determined by one- and two- period rates vectors 
R = [Ri-i^i]i and i?*^^) = [Ri-2,i\i, the system is closed. We can write a vector SDE 
which is a vector diffusion for all the one- and two- period rates under any of the Ci_i,i 
measures: 



d 



R 



diag(/x(i?,i?(2)),0(i?, R^^^)) 



R 



dt + diag((T, 



R 



d 



Z' 



At this point a Monte Carlo simulation of the process, based on a discretization 
scheme for the above vector SDE is possible. One only needs to know the initial CDS 
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rates R{0), i?^^^(0), which if not directly available one can build by suitably stripping spot 
CDS rates. Given the volatilities and correlations, one can easily simulate the scheme 
by means of standard Gaussian shocks. 

If C is the Gholesky decomposition of the correlation Q {Q = CC with "C" lower 
triangular matrix) and 1^ is a standard Brownian motion under Cj_i,j, we can write 



d 



R 



= diag(//(i?,i?(^))>(i?,i?(2))) 



R 



dt + diag((T, v) 



R 



C dW (11) 



The log process can be easily simulated with a Milstein scheme. 



5 An approximated model with a single family of 
rates and proof of the main result 

Consider now the approximation 



Rj{t) '.— Lgd 



E[D{t,T^)l{T,.,<r<T,}\J't. 



= Lc 



E[D(t,T,)l{.>r,_,}|^f]-E[D(t,T,-)l{,>r,}|.^t] 
ajQir > t\J^t)P{t,Tj) 



E[L>(t,r,_i)lu>r,_,} 


Dit,Tj)/Dit,Tj_^) 


JS]-E[L>(i,7;-)l{.>r,}|^t] 


ajQ{T > t\Tt)P{t:Tj) 



= Lc 



At this point we approximate the boxed ratio of stochastic discount factors with the 
related zero-coupon bonds, obtaining 



jGD 



E[D(i,r,_oi{.>T,_,}l^t] 


P{t,T,)/Pit,T,^^) 


-E[D(i,r,)l{,>r,}|^d 


a 





= Lc 



P{t, Tj_,)P{t, Tj)/P{t, T,_i) - P{t, Tj) Lgd 



ajP{t, Tj) 



P{t,T^-i) 



P{t,T, 



-^l]^R^{t) 



- 1 



a,- \{l + a,F,mP{t,T,] 

where F is the forward LIBOR rate between Tj_i and Tj. This last definition can be 
inverted so as to have 



P{t,Tj) VLgd^-^ ' -J^- ' -.-.V"//' - (12) 

as long as i? > 0, provided that Fj{Q) > as should be. This means that we are free 
to select any martingale dynamics for Rj under Q^"^'-', as long as Rj remains positive. 
Choose than such a family of R as building blocks 



dRi{t) = ai{t)Ri{t)dZl{t), for all i 
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and define the P by using ( fT2|) to obtain inductively P(t,Tj) from P(t,T,_i) and 
from Rj. This way, the numeraires P become functions only of the i?'s, so that 
now the system is closed and all one has to model is the one-period rates R 
vector. No need to model two-period rates in this framework. 

In this context the change of numeraire becomes 



dZ^ = dZ' - pDC 



Hfe 



dt = dZ' - pDC In 



Oh 
Lgd 



dt 



H (f^^i?, + l)(l + «,F,(0)) 
\h=j+i 



dZ'- pDCln 

i 

dZ^-pY, DCln((^^, + l)(l + a;,F,(0)))' 



dt 



dt 



h=j+i 



dZ'-p y DCln f f P^Rh + l)) dt = dZ'-p V . — BCiRh)' 

h=j+l + — 



dt 



so that we can write 



dZi = dZl-Y. PK^^^^dt 
h=j+l + — 

from which we have the dynamics of Ri under Q-': 

dRi = a.RidZi = (TiR, I dZi + J2 P^'^^^^^^^U^dt I =: R,{fxliR)dt + aidZl) (13) 

Consider the drift term in the last formula. If we compute -E-'~^'-' [i?j(Tj„i)] we obtain 

E^-^^^[R,{T,_^)] ^ R,iO)expy^' 'jll{R{0))du 

= ^,(0)exp| V .^—ML^p f' \,(u)ak{u)du\{U) 
Uf^i/?fc(0) + LGD/afc Jo J 

and, if we take volatilities a to be constant, we have 



\k=j+l 



Under independence between intensities and interest rates (and in particular under 
deterministic intensities, which are a common assumption when stripping one-period 
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CDS rates from multi-period ones), by definition of Rj it is easy to show that at time 0, 
both the original Rj{0) and the approximated Rj{0) are given in terms of the survival 
probabilities as 

i?,(0) = ^,(0) = LcB/a, [^^^ - l) (15) 
and hence (flU) is reduced to 



(16) 

It follows that the convexity effect vanishes if the ratio is close to one. 

Now, based on the approximated dynamics ( fT3ll and the related expectation above, 
we prove the main result of the paper, i.e. Proposition 12. 4[ 

Proof. To prove the proposition, we compute the price of the premium leg as 

6 

J2 ajEo[D{0,T^)l{,^T,}Rj-i,MTj-i)] = ■■■ 

j=a+l 

The first approximation we consider is applied to Rj^ij+ciTj^i), so that 
Then by substituting this in the premium leg expression we have 

b j+c 

E E"^-^'(o)^o[^(0'^^)w.}^^(^.-i)] = 

j=a+l i=j 
b j+c 

= E Y.''Mmo[Dio,TJ)R^iT,.,)E{^,^T,}\J'T,)] 

j=a+l i=j 



j=a+l i=j 

= E J2''mc,-iAmi'''mT,-i)]= E E«,^^(o)p(o,r.)Er"[^.(^.-i)] = 

j=a+l i=j j=a+l i=j 

where we have applied the change of numeraire, moving from the risk neutral numeraire 
B to the numeraires Cj^i/s. The last expected value can be computed based on (fT6!) . 
By substituting the expected value expression, we obtain the final formula. 

□ 
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6 A few numerical examples 

We report input data and outputs for a name with relatively large CDS forward rates. 
We consider the FIAT car company CDS market quotes as of December 20, 2004. Since 
in Brigo and Alfonsi (2003) we have some evidence on the fact that CDS prices depend 
very little on the correlation between interest rates and credit spreads, when stripping 
credit spreads from CDS data we may assume independence between interest rates and 
credit spreads. This leads to a model where it is easy to strip default probabilities 
(hazard rates) from CDS prices, as hinted at again in Brigo and Alfonsi (2003). Using 
this independence assumption, we strip default (or survival) probabilities from CDS 
quotes with increasing maturities. 

6.1 Inputs 

We take as inputs the following Fiat CDS rates and use mid quotes 



n 




pASK 


lY 


99.9 


175.57 


2Y 


172.5 


231.38 


3Y 


243.73 


286.13 


5Y 


348.85 


366.54 


7Y 


380 


410 


lOY 


395.16 


412.73 



We take Rec = 0.4 (so that Lgd = 0.6). The input zero coupon curve, and the 
survival risk-neutral probabilities stripped from the above CDS quotes are reported in 
Appendix 1 below. 

6.2 Outputs 

We start by giving a table for 

Conv(cr, p) := CDScMa^b^c(0, Lgd, a, p) — CDScMa^{,^c(0, Lgd; p = 0). 

The first term is computed by assuming the volatilities o"j of forward one-period CDS 
rates Ri to have a common value a and the pairwise correlations pij to have a common 
value p. This first term is then given by formula ([T]). The second term is the simpler 
value ([2]) where no correction due to CDS forward rate dynamics is accounted for. This 
difference then gives us the impact of volatilities and correlations of CDS rates on the 
CMCDS price. The difference is always positive, similarly to what happens to analogous 
constant maturity swaps in default free markets under similar conditions on volatilities 
and correlation. It is the impact of "convexity" on the CMCDS valuation. We take 
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a = 0, b — 20 (5y final maturity) and c = 21 (which means we are considering non- 
standard 5y6m CDS rates in the CMCDS premium leg, c = 19 would amount to a 5y 
CDS rate). 
We obtain 



Conv((T, p) 


p: 0.7 0.8 0.9 0.99 


a: 0.1 
0.2 
0.4 
0.6 


0.000659 0.000754 0.000848 0.000933 
0.002662 0.003047 0.003435 0.003784 

0.011066 0.012742 0.014442 0.015995 
0.026619 0.030964 0.035464 0.039652 



The "convexity difference" increases with respect both to correlation and volatihty, 
as expected. 

The next table reports the so called "participation rate" 0a,fe,c(o", p) for a CMCDS 
with final Th = by {a = 0,b = 20, recalling that resets occur quarterly), with 5y6m 
constant maturity CDS rates (c = 21), 



00,20,21 (o-, p) 



E-li«,P(0,T,)i?o,2o(0) 



"premium leg CDS" 
"premium leg CMCDS" = E?=i ^'^^(0, T,)^^>t,}Rj-i,MT,-i)] ' 



The CMCDS premium leg is computed with our approximated market model based 
on one-period rates R. As we see from the outputs, the participation rate increases with 
volatility and correlation, as is expected from the "convexity adjustment" effect. 



00,20,21 (O-, P) 


p: 0.7 0.8 0.9 0.99 


a: 0.1 
0.2 
0.4 
0.6 


0.71358 0.71325 0.71292 0.71262 
0.70664 0.70532 0.704 0.70281 
0.67894 0.67368 0.66842 0.66368 
0.63302 0.62128 0.60957 0.59907 



Finally, we fix volatilities and correlations and check how the patterns change when 
changing final maturity Tb = Ti. We consider the following quantities at time and with 



"Constant maturity rate" i?j_i^j+c(0) . ^ ^ 
* "standard rate" -Ro,fe(0) ' ' ' 

_ E^-^'-[L>(O,T0l{.>r.}i?.-M+cm-i)] 

p{o,T,)RoM ' ' 

Er^'^[i^(o,r,)i{.>Toi?»-i,^+c(ri-i)] .^^ ^ 
p(o,r,)i?i_i,,+c(o) ' ' 

"premium leg CDS" _ Ej=i Q:j-P(0, ^j)-Ro,»(0) ^ 
~ "premium leg CMCDS" " «,-P(0, r,)i?,_i,,+,(0) ' ' ~ ' ' ' ' ' 
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"premium leg CDS" Ej=i «j-P(0, ^j)-Ro,»(0) 

~ "premium leg CMCDS with convexity" ~ jy.^^ q;X-''^[L>(0, Tj)l{r>T,}Rj-i,j+c{Tj-i)] ' 





Xi 


Vi 








(pi 


















a 


= 0.4; 
















P 


= 0.9; 


1 


0668 


1.0668 


1 





37773 


0.37773 






1 


1288 


1.1359 


1.0063 





36281 


0.36162 


Rec 


= 0.4; 


1 


1914 


1.2075 


1.0135 





35281 


0.35039 




a=0; 


1 


2525 


1.2792 


1.0214 





34359 


0.33993 


c 


= 20; 


1 


3107 


1.3495 


1.0297 





33512 


0.33024 


b 


= 20; 


1 


3673 


1.4193 


1.038 





34187 


0.33548 






1 


4171 


1.4826 


1.0462 





36905 


0.36064 






1 


4515 


1.53 


1.0541 





40755 


0.39664 






1 


4716 


1.5622 


1.0616 





45262 


0.43881 






1 


4798 


1.5818 


1.0689 





49477 


0.47785 






1 


4837 


1.5979 


1.0769 





52661 


0.50671 






1 


4905 


1.6175 


1.0852 





55072 


0.52799 






1 

i 


A nnn 

4yyy 


i.D4Uo 


i.UyoD 





obyol 


f\ ^A OOA 






1 


5122 


1.666 


1.1018 





58674 


0.55846 






1 


5236 


1.69 


1.1092 





60704 


0.57574 






1 


5275 


1.706 


1.1168 





62715 


0.5928 






1 


5274 


1.7174 


1.1244 





64681 


0.60938 






1 


5249 


1.7236 


1.1303 





67017 


0.62939 






1 


5106 


1.7173 


1.1368 





69254 


0.64843 






1 


4924 


1.7047 


1.1422 





71589 


0.66842 







The Xj's measure how the constant maturity CDS rate differs multiplicatively from 
the standard CDS rate, so they are a measure of how the constant maturity CDS differs 
from a standard CDS in the premium rate paid at each period. We find an increasing 
pattern in Tj as partly expected from the fact that the input CDS rates R^l^'^^^ are 
increasing with respect to maturity Ti,. 

The yiS measure the same effect while taking into account "convexity" , i.e. future 
randomness of the payoff and correlation. The t/^'s would reduce to the x^'s if correlations 
p were taken equal to 0. The y maintain the increasing pattern with respect to Tj. 

The Zi^s measure the multiplicative impact of "convexity", in that they are due to 
contributions stemming from volatilities a and correlations p of CDS rates. The impact 
is increasing with maturity Tj, as expected from the sign in the exponent of the convexity 
adjustments and from the positive signs of correlations (and volatilities). 

Finally, as seen above, the V'i's are the so called "participation rates" for different 
terminal maturities T^. They give the ratio between the premium leg in a standard CDS 
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protecting in [0, Tj] and the premium leg in CMCDS for the same protection interval when 
ignoring the convexity adjustment due to correlation and volatilities. The 0j's are the 
participation rates computed when taking into account convexity due to volatilities and 
correlations. We have seen a particular participation rate (j) earlier. In the table above 
for we obtain an initially decreasing pattern followed by a longer increasing pattern 
for both ip and 0. Notice that, on the longest participation rate, in the last row of the 
related table, convexity has an impact moving from a 71.59% participation rate when not 
including "convexity" (ignoring correlations and volatilities) to a 66.84% participation 
rate when including convexity. There is a 4.8% difference in the participation rate of 
this FIAT 5y-5y3m CMCDS with correlations set at 0.9 and volatilities at 40%. 



7 Further work 

In further research we need to propose a realistic dynamics for one- and two- period rates 
that completely specifies the market model, along the guidelines given in this paper and 
in Brigo (2004). Then we may test the market formula proposed here against Monte 
Carlo simulation of the exact dynamics. Moreover, examining the formula outputs for 
p matrices with more realistic decorrelation patterns and for different names can be 
appropriated. 
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Appendix 1. Input data and survival probabilities 

The input zero coupon curve, and the survival risk-neutral probabilities stripped from 
FIAT CDS quotes are: 
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24444 





25556 





25556 





25278 




0.25 





25556 





25556 





25278 




0.25 





25556 





25556 





25278 





25278 





25556 





25556 





25278 




0.25 





25556 


n 
u 







25278 




0.25 





25556 





25556 



T- 








0.24444 




0.5 


0.75556 


1 


0083 


1 


2583 


1 


5139 


1 


7694 


2 


0222 


2 


2722 


2 


5278 


2 


7833 


3 


0361 


3 


2889 


3 


5444 




3 8 

O.O 


4 


0528 


4 


3028 


4 


5583 


4 


8139 


5 


0667 


5 


3167 


5 


5722 


5 


8278 



0.99994 
0.99459 
0.989 
0.98309 
0.97709 
0.97098 
0.96458 
0.958 
0.95133 
0.94448 

0.9373 
0.92995 
0.92251 
0.91502 
0.90731 
0.89929 
0.89139 
0.88373 
0.87544 

0.8673 
0.85906 
0.85085 
0.84255 
0.83417 



Q(r > T,) 
0.99994 
0.99429 
0.98856 
0.98279 
0.97712 
0.97155 
0.96433 
0.95409 
0.94108 
0.92552 

0.9086 
0.89227 
0.87669 
0.86165 
0.84618 
0.82931 
0.81203 
0.79449 
0.77495 
0.75531 
0.73503 

0.7142 
0.69403 
0.67559 
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D. Brigo, constant maturity CDS pricing with the Market model 





Ti 


P(O.T,) 


Q(r > Ti) 


0.25278 


6.0806 


0.82572 


0.65879 


0.25 


6.3306 


0.81735 


0.64351 


0.25556 


6.5861 


0.80892 


0.62968 


0.25556 


6.8417 


0.80035 


0.61709 


0.25278 


7.0944 


0.79182 


0.60594 


0.25278 


7.3472 


0.78344 


0.59601 


0.25556 


7.6028 


0.77494 


0.58651 


0.25556 


7.8583 


0.76641 


0.57685 


0.25278 


8.1111 


0.75794 


0.56715 


0.25 


8.3611 


0.74977 


0.55744 


0.25556 


8.6167 


0.74141 


0.5474 


0.25556 


8.8722 


0.73303 


0.53726 


0.25278 


9.125 


0.72474 


0.52713 


0.25 


9.375 


0.7168 


0.51705 


0.25556 


9.6306 


0.70869 


0.50667 


0.25556 


9.8861 


0.70041 


0.49601 


0.25278 


10.139 


0.69241 


0.48565 


0.25 


10.389 


0.6849 


0.4756 



